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Abstract 

Quasi-particle model of quark gluon plasma is the statistical mechanics of particles with medium de- 
pendent mass, related to plasma frequency, which was proposed to describe the thermodynamics of the 
medium itself. At relativistic limit the plasma frequency depends on number density and temperature. 
The number density is a thermodynamic quantity of the medium which in turn depends on plasma fre- 
quency. Hence, one need to solve this problem self-consistently, instead of using perturbative expressions 
for plasma frequency. Here we carry out such a self-consistent calculations using our, recently developed, 
a new formulations of quasiparticle model. By adjusting a single parameter for each system, a remarkably 
good fit to lattice QCD results are obtained for 2, 3 and (2-1-1) flavor quark gluon plasma systems, first, 
with zero chemical potential. Then, it is extended to systems with finite chemical potential and fits very 
well the lattice results without any new parameter. 
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1 Introduction : 

Recently proposed self-consistent quasiparticle model for quark gluon plasma (qQGP) |T] is extended to 
2-flavor (2-f), 3-flavor (3-f) and (2+l)-flavor ((2+l)-f) QGP and is found to explain very well the lattice 
gauge theory (LGT) QCD results 2 . Our model differs from other qQGP models [3] in many aspects. 
Firstly, our formulation is based on the standard statistical mechanics (SM) [4J, without the need of any 
reformulation of SM [S]. Here we start from energy density (e) and particle number density (n), well defined 
in grand canonical ensemble (GCE) and then all other thermodynamic (TD) quantities, including pressure, 
are derived from e, n and using the standard TD relations. Hence, there is no TD inconsistency in our model 
[6]. Secondly, the thermal mass of quasiparticle, in our model, is assumed to be function of density which 
in term depends on thermal mass and hence one need to solve it self-consistently [T] with the requirements 
that at high temperature thermal masses go to that of QCD perturbative results. Note that all other qQGP 
models [3j used QCD perturbative thermal masses which is appropriate to ideal system, because one uses 
ideal thermal propagators in these calculations, but the system is non-ideal. Only at high temperature it 
may be approximated to ideal system. Thirdly, there is only single system dependent adjustable parameter 
and gives a surprisingly good fit to LGT results. We will discuss and compare 3 models, in our formulation, 
based on 3 forms of thermal mass. Model-I is based on QCD perturbative thermal mass without self- 
consistent calculation. Model-H is based on our self-consistent calculation, using the simpler, approximate 
dispersion relation {ujk Vs k). In model-HI, we repeat the self-consistent calculations of Model-H, but using 
the dispersion relation of hard-thermal-loop (HTL) resumed QCD perturbation theory. We will see that the 
results of Model-I, is not satisfactory, but results improve considerably as we go from Model-I to Model-HI, 
because of our self-consistent calculations. 

2 New formulation of qQGP: 

As in any qQGP model, here also, we assume that the thermal properties of interacting real particles are 
modeled by non-interacting quasiparticles with additional thermal mass which depends on thermodynamic 
quantities due to collective properties of the medium. Following the the standard SM, in GCE, e and n 
are the only TD quantities, defined as the average of the energy and the number of particles. All other TD 
quantities, including pressure, need to be derived [71 S]. In usual SM, the masses of particles are zero or 
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constant and hence one can derive the expression for the pressure from e, using TD relations, and find that 
the pressure is equal to the logarithm of the partition function [7]. But, here in QGP it is not true, since 
the thermal mass depends on TD quantities which in turn depends on thermal mass through a functional 
relationship and hence whole problem need to be solved self-consistently. 
The energy density is [1] 



9f 



dk 



(1) 



where g/ is the degeneracy and =F refers to bosons and fermions. z is the fugacity. Note that, here we 
have assumed that the whole thermal energy is used to excite quasiparticles and hence the vacuum energy 
contribution is neglected. A detailed discussion on the vacuum energy contribution and TD consistency 
relation was commented in Ref. [6]. Similarly, the number density is 



9f 



dk k^ 



1 
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Single particle energy or tOk depends on thermal masses and momentum k. The exact expression for ivk may 
be obtained by solving the approximate dispersion relation, obtained using HTL calculations [i8j|. 
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for quarks, where LUp and ruf are plasma frequencies related to thermal masses. Both of them are, in general, 
functions of temperature (T) and chemical potential (/i). First, we consider QGP system with /i = and at 
the end with /i 7^ 0. Above dispersion relation may be further approximated to a simpler form. 



and 

for gluons and quarks respectively, valid at high temperatures. The thermal masses are defined as 
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TO^ = -w.^ and = 2m ^ , 
respectively, for massless particles. For massive quark with mass mo, m^ may be modified as [5] 

m^ = (mg + TO/)^ + mj . 
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Figure 1: Plots of PjT'^ as a function of T/Tc from the Model-I and lattice results [2] (symbols) for pure 
gauge, 2-f, 3-f and (2+l)-f QGP. 



First, let us consider our Model-I where, just like in other qQGP models, we take 

(6 + n/) , 



18 



and 



9^ 
6 



(9) 



(10) 



where g is the QCD coupling constant and n/ is the number of flavors. These are perturbative results, valid 
for very high temperature, but are modified to take account of non-perturbative effects by using appropriate 
running coupling constant, as{T) = j^. Different authors [3j use different phenomenological models of as 
and even thermal masses, with few adjustable parameters. Here, in our Model-I, we use 2-loop approximate 
running coupling constant, 



3(153 - 19n/) ln(2 ln(r/ At)) \ 
(33-2n/)2 ln(T/AT) ) ' 



(11) 



(33 - 2n/)ln(r/AT 

where At is the only one adjustable parameter in our model. Results are presented in Fig. 1 and Fig. 2 
The general expression for the energy density is 
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Figure 2: Plots of e/T"* as a function of T/Tc from the Model-I and lattice results [2] (symbols) for pure 
gauge, 2-f, 3-f and (2+l)-f QGP. 

where gg — 16 is the degeneracy of gluons and n^^^ , 'nf'Jf are the effective number of flavors corresponding 
to light quarks [u and d) and strange quark respectively, n'j^^ ^ fo'' rnassless quarks and is less than n f 
for quarks with finite mass |2 . 

The pressure is obtained from the thermodynamic relation [T] 



dP 

e = T- P 

dT 



(13) 



on integration with one integration constant which may be chosen such that P Pg, the Stefan-Boltzmann 
limit, as r — > oo. However, it is difficult to adopt it in our numerical work and hence we fix it at T = Tc to 
LGT results. 

3 Self-consistent qQGP models - Model-II and Model-Ill: 

To improve Model-I, in Model-II, we replace the perturbative expression for plasma frequencies, Eq. Q and 
Eq. (jlOp . by density dependent expression, 
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for gluons and 
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which is motivated from similar work in rclativistic electron-positron plasma [21 llOj where we know that the 
plasma frequency, at relativistic limit, is proportional to n/T. We fix the constant of proportionality a^, aq 
and hq by demanding that as T — > oo, ujp and nif goes to the corresponding perturbative results. Since now 
thermal masses depends on density, the expression for density becomes a functional, given by 
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where x = k/T and uj = ujk/T which depends on x, ^ and g^{T). 
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Figure 3: Plots of P/T^ as a function of T/Tc from the Model-II and lattice results [2] (symbols) for pure 
gauge, 2-f, 3-f and (2-Hl)-f QGP. 



So we need to solve density equation self-consistently to get the density and then we may calculate plasma 
frequencies which may be used to evaluate e, P and other TD quantities using Eq. and Eq. It 
may be convenient to work with 

/•oo -1 

(17) 



for gluons, 



for quarks and 
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Figure 4: Plots of e/T'^ as a function of T/T^ from the Model-II and lattice results [2] (symbols) for pure 
gauge, 2-f, 3-f and (2+l)-f QGP. 

for antiquarks, where fg = ^ fq = ■j^ ^ and /| = ^ = 6n|^''^''^ is the degeneracy of 

quarks and antiquarks. It is easy to see that at high temperature ^ and hence ^ 2 C(3) and 
fq=fq^ 2 77(3) for /i = systems, appropriate to ideal plasma. This high temperature limit is taken to 
evaluate the coefficients a^, Og and bq by comparing with perturbative results for plasma frequencies. Results 
are presented in Fig. 3 and Fig. 4. 

Next we further improve Model-II by using the full HTL approximate dispersion relation, Eq. ([3]) and 
Eq. Q, instead of the approximate dispersion relation. First, one need to solve HTL dispersion relation for 
LOk and then carry out the self-consistent calculations as done in Model-II. Results are presented in Fig. 5 
and Fig. 6. We see that there is only a marginal improvement, especially near to T < 2Tc, as we go from 
Model-II to Model-Ill. Hence the dispersion relation of Model-II is a good approximation to that of HTL 
calculations. 

4 Self-consistent qQGP with /x 7^ 0: 

So far we have discussed QGP systems with /i = and found to explain LGT results very well using our 
self-consistent model. Let us apply our model to the system with ji ^ {). For ji ^ {), thermal masses and 
running coupling constant need to be modified as a function of T and jj,. In our self-consistent model we 
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Figure 5: Plots of P/T'^ as a function of T/Tc from the Model-Ill and lattice results [2] (symbols) for pure 
gauge, 2-f, 3-f and (2+l)-f QGP. 



demand that as T oo, mj: goes to QCD perturbative (T, /i) dependent functions, 

,2 



2 g'T' , r 

"If = — (1 



The running coupling constant may be modified by replacing T/ At in Eq. pip by 
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(21) 



using the results of Schneider [TT . Similar /i dependent modifications also used by Letessier, Rafelski [12], 
but instead of (1.91/2.91)^ factor they used I/tt^ which also gives similar results except near Tc- First, we 
solve the density equations, Eq. (jlSp and Eq. (|19p . for quarks and antiquarks self-consistently and then 
subtract each other to get the net quark density Uq/T^ . Note that now the asymptotic values of fq and /g, 
needed to evaluate m/, are not 2?7(3), but given by 
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Figure 6: Plots of e/T"^ as a function of T/Tc from the Model-Ill and lattice results [2] (symbols) for pure 
gauge, 2-f, 3-f and (2+l)-f QGP. 

and both the expressions approach 277(8) for z = 1. Results of our model, using Model-II, is plotted in Fig. 
7 for riq/T^ and fits LGT results without the need of any new parameters and Kt is fixed earlier for /i = 
system. 

5 Results: 

The analysis using Model-I is straight forward and the parameter to = is adjusted to get the best fit 
to LGT results on pressure. Using this ta we evaluate other TD quantities like e. The fitted values of io 
are 0.35, 0.4, 0.2, 0.2 for pure gauge, 2-f, 3-f and (2+l)-f QGP respectively and the results are plotted in 
Fig.l and Fig. 2. In the case of massive quarks, effective number of quarks flavors, n^/^ , is less than n/ and 
may be obtained from the Ref. :2] as n^Jf = 1.8344 and nlff = 0.8275. The LGT data ^2: for the systems 
with finite flavor were multiplied by a factor 1.1 to take into account of the uncertainties in the continuum 
extrapolation. The fit to LGT results are not satisfactory for T < 2.5 Tc, but fits well at higher temperature 
as it should. Note that the same model was also used by other qQGP studies [3] and obtained very good fit 
to LGT results, but adjusting 3 or more parameters, where as here we have only single adjustable parameter. 
In the case of Model-II, we solve the equation for density, self-consistently and use it to evaluate the pressure. 
Again, the parameter to is adjusted to fit LGT results on pressure and then the energy density is evaluated. 
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Figure 7: Plots of riq/T^ as a function of T/T^ from the Model-II and lattice results ^ (symbols) for 2-f 
QGP. 

Now the fitted values of are higher and are 0.7, 0.6, 0.35, and 0.45 for pure gauge, 2-f, 3-f and (2+l)-f 
QGP respectively and the fit to LGT results. Fig. 3 and Fig. 4, improves a lot compared to Model-I. Similar 
self-consistent calculations is used in Model-Ill using the HTL dispersion relation, and the results are plotted 
in Fig. 5 and Fig. 6. The fitted values of to are 0.74, 0.68, 0.4 and 0.55 for pure gauge, 2-f, 3-f and (2+l)-f 
QGP respectively. We see further improvement in the results, especially near T < 2Tc, and the values of to 
are slightly higher. These values of to may give better values of At, since Model-Ill is more accurate than 
Model-I and Model-II. In Fig. 7, results of our Model-II for 2-f QGP system at finite /i is compared with 
LGT results for quark density and reasonably good fit is obtained without any new parameter. Of course, 
very close to T = Tc and for T < Tc, our quasiparticle model, which is based on plasma collective effects, 
may not explain the LGT results. One may need to include the effects of vacuum energy or effects of hadron 
resonance gas as done in Ref. [13]. 

6 Conclusions: 

Here we have studied the thermodynamics of pure gauge, 2-f, 3-f and (2+l)-f QGP systems with our newly 
developed, self-consistent quasiparticle model using HTL approximate dispersion relation and it's approxi- 
mate simpler dispersion relation. By adjusting single system dependent adjustable parameter, to = At/T^, 
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we obtained a remarkable good fits to LGT results, without the need of too many adjustable parameters 
as required by other qQGP models [3 . For comparison, we have also studied these systems without using 
the self-consistent calculations, but fails to fit LGT results for T < 2.5 Tc- From our calculation on /i 7^ 
system, if LGT results are true, we need further more accurate dispersion relation than HTL approximate 
results to improve our model. It may be interesting to adopt this idea of self-consistent calculations in finite 
temperature QCD and solve the QCD at finite temperature without the need of qQGP. 
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